NASA/CR-1999-209547 



Modeling of Longitudinal Unsteady 
Aerodynamics of a Wing-Tail Combination 

Vladislav Klein 

The George Washington University 

Joint Institute for Advancement of Flight Sciences (JIAFS) 

Langley Research Center, Hampton, Virginia 


September 1999 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated 
to the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key 
part in helping NASA maintain this 
important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the lead center for 
NASA's scientific and technical information. 
The NASA STI Program Office provides 
access to the NASA STI Database, the 
largest collection of aeronautical and space 
science STI in the world. The Program Office 
is also NASA's institutional mechanism for 
disseminating the results of its research and 
development activities. These results are 
published by NASA in the NASA STI Report 
Series, which includes the following report 
types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results 
of NASA programs and include extensive 
data or theoretical analysis. Includes 
compilations of significant scientific and 
technical data and information deemed 

to be of continuing reference value. NASA 
counterpart or peer-reviewed formal 
professional papers, but having less 
stringent limitations on manuscript 
length and extent of graphic 
presentations. 

• TECHNICAL MEMORANDUM. 

Scientific and technical findings that are 
preliminary or of specialized interest, 
e.g., quick release reports, working 
papers, and bibliographies that contain 
minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. 

Collected papers from scientific and 
technical conferences, symposia, 
seminars, or other meetings sponsored or 
co-sponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to 
NASA's mission. 

Specialized services that complement the 
STI Program Office's diverse offerings include 
creating custom thesauri, building customized 
databases, organizing and publishing 
research results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home 
Page at http:Hwzvw.stiMasa.gov 

• Email your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA STI 
Help Desk at (301) 621-0134 

• Telephone the NASA STI Help Desk at 
(301)621-0390 

• Write to: 

NASA STI Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076-1320 


NASA/CR-1999-209547 



Modeling of Longitudinal Unsteady 
Aerodynamics of a Wing-Tail Combination 

Vladislav Klein 

The George Washington University 

Joint Institute for Advancement of Flight Sciences (JIAFS) 

Langley Research Center ; Hampton , Virginia 


National Aeronautics and 
Space Administration 


Langley Research Center 
Hampton, Virginia 23681-2199 


Prepared for Langley Research Center 
under Cooperative Agreement NCC1-29 


Available from: 


NASA Center for AeroSpace Information (CASI) 
7121 Standard Drive 
Hanover, MD 21076-1320 
(301) 621-0390 


National Technical Information Service (NTIS) 
5285 Port Royal Road 
Springfield, VA 22161-2171 
(703) 605-6000 


Summary 


Aerodynamic equations for the longitudinal motion of an aircraft with a horizontal tail were devel- 
oped. In this development, emphasis was given on obtaining model structure suitable for parameter esti- 
mation from experimental data. The resulting aerodynamic models included unsteady effects in the 
form of linear indicial functions. These functions represented responses in the lift on the wing and tail 
alone, and interference between those two lifting surfaces. The effect of the wing on the tail was formu- 
lated for two different expressions concerning the down wash angle at the tail. The first expression used 
the Cowley-Glauert approximation known as “lag-in-downwash,” the second took into account growth 
of the wing circulation and delay in the development of the lift on the tail. Both approaches were dem- 
onstrated in two examples using the geometry of a fighter aircraft and a large transport. It was shown 
that the differences in the two downwash formulations would increase for an aircraft with long tail arm 
performing low-speed, rapid maneuvers. 
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Introduction 

Recently there has been an increased interest in identification of aircraft aerodynamic models that 
include unsteady aerodynamic terms. The majority of identified models were related to the longitudinal 
motion of tailless aircraft, see reference 1 to 3, where the aerodynamic models had relatively simple 
form. In formulating aerodynamic model equations for a wing-tail combination, however, more compli- 
cated model forms can be expected. This problem was first addressed in the early twenties by Cowley 
and Glauert (ref. 4). They realized that there was a time lag before the wing disturbance reaches the tail. 
They assumed that the downwash associated with a change in lift is equal to the corresponding steady 
value but the effect at the tail is delayed by the time for the airplane to travel a distance equal to the tail 
arm. The change of the lift on the tail due to downwash was included in the damping-in-pitch deriva- 
tive. Now this change is interpreted as a tail contribution to the acceleration derivatives, i.e., the deriva- 
tive of the lift and pitching moment with respect to rate of change in the angle of attack. 

The investigation of the downwash angle was extended by R. T. Jones and Fehlner in reference 5 by 
considering both the growth of wing circulation and the delay in the development of lift by the tail. 
They presented an expression for the downwash indicial function associated with a change of the lift on 
the wing. Considerable effort to the rigorous explanation of the downwash both qualitatively and ana- 
lytically was given by Tobak in reference 6. He developed indicial functions for the lift and pitching 
moment of a wing-tail combination and presented numerical results for several representative cases in 
supersonic flight regimes. The theoretical results were then applied to the longitudinal stability analysis 
of an aircraft. 

In reference 7 a theoretical analysis of the downwash due to changes in the lift of the wing was pre- 
sented. In this analysis the wing was approximated either by a single horseshoe vortex or by a set of 
these vortices. The result showed that for specific changes in the lift the downwash would depend on the 
tail arm, wing aspect ratio and its span. Similar to reference 7, a simple vortex system was used in refer- 
ence 8 to model unsteady aerodynamic effects into the longitudinal equations of motion. It was found 
that the solution of resulting integro-differential equations led to lengthy computations. The calculations 
also showed that the unsteady aerodynamics could be adequately accounted for when the unsteady 
effects were included only in the downwash development. This finding permitted reformulating the 
equations so that they appeared suitable for parameter estimation using flight test data. 

This report is an extension of previous research in references 5, 7, and 8. Its purpose is to develop 
aerodynamic equations for the longitudinal motion of an aircraft with horizontal tail. These equations 
include unsteady effects modeled by linear indicial functions. The emphasis is given to obtaining model 
structure suitable for model identification using experimental data. It means that the postulated model 
should explain the aerodynamics of wing-tail combination with sufficient accuracy and, at the same 
time, should be simple enough to provide good conditions for parameter identifiability. The report starts 
with the formulation of equations for the lift and pitching-moment coefficient and introduction of indi- 
cial functions in these equations. Then the development of an indicial function relating the change in the 
lift of a wing due to downwash is presented. The downwash angle during the harmonic motion of an air- 
craft is compared with that based on the traditional “lag-in-downwash” factor contained in the accelera- 
tion stability derivatives. The variation of the downwash angle with the geometry and frequency of the 
harmonic motion is also briefly discussed. Next, the final form of aerodynamic model equations is pre- 
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sented. Numerical examples using the geometry of a fighter and transport aircraft demonstrate possible 
forms of indicial functions using both the ‘iag-in-downwash” formulation and that with unsteady terms. 
The report is completed by concluding remarks. 

Indicial Functions 


In the following development the lift and pitching moment will be considered as functions of the 
angle of attack and pitching velocity, i.e., 

C a = C a (a, q), a-Loxm 

Then the aerodynamic model equations can be formulated as 


t t 

C a U) = C a (Q) + jc a Jt -'C)Ct(T)dT + (t-T)q(x)dx 

0 0 


( 1 ) 


Where C (0) is the value of the coefficient at initial steady-state conditions, and C (t) and C (/) 

a a a q 

are the indicial functions defining the responses in C to unit step in a and q respectively. Replacing 
the second indicial function in ( I ) by its steady value C (oo) , equation ( 1 ) is simplified as 

t 


c a (t) = C a (0) + J^c (oo)qit) + jc a Jit-x)d(x)dx 

0 

When the deficiency functions 


( 2 ) 



U) = c a («>)-c a (0 

a a 


are substituted into (2) the equations for the aerodynamic coefficients take the form 


C a (t) = C fl (0) + C flo (oo)a(0 + ^C (oo)^(/)-jF flo (/-T)a(T)rfr 

o 


(3) 


(4) 


where C (•») and C (°°) are the rates of change with a and q of the coefficients C , and C . 

11 a U q L 


The indicial functions C Q ( t ) include combined responses of the wing and tail, and interference 
effects between those two lifting surfaces. Because of linear aerodynamics the resulting indicial func- 
tions are given as a sum of their components. For the indicial function C L (t) these components will be 
as follows: 


( 1 ) the response in the lift of the wing to a unit step in angle of attack of the wing, while the angle of 
attack of the tail remains at zero; 

(2) the response in the lift of the tail to a unit step in angle of attack of the tail, while the angle of 
attack of the wing remains at zero; 
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(3) the response in the lift of the tail to a unit step in the angle of attack of the wing, while the angle 
of attack of the tail remains at zero; 


(4) the responses in the lift of the wing to a unit step in the angle of attack of the tail, while the angle 
of attack of the wing remains at zero. 

The first two components are represented by the response of an isolated wing and isolated tail 
respectively. As pointed out in reference 6, these representations are true for supersonic speeds. For 
subsonic speeds, however, these are good approximations. The third component expresses the lift on the 
tail due to a change in the downwash induced by the lift of the wing. It is usually a significant contribu- 
tion to the resulting pitching moment of an aircraft. Finally, the last component is zero for supersonic 
speeds. For subsonic speeds it can be neglected. A summary of all four components, their conditions 
and interpretations is given in table I . 


The expressions for the indicial function C Lw (t) of a wing alone were given by several authors. 
In reference 5, the indicial function was formed as°a sum of exponential 


C 


L\\\ 


( t ) = 


c* 0 + 




V 


(5) 


where the values of parameters c 0 , Cj and X . depend on the wing geometry. For a wing with small 
aspect ration equation (5) was further simplified as 




( 6 ) 


Equation (6) was used in references 1 and 3 for aerodynamic model identification of tailless aircraft 
from wind tunnel, small amplitude oscillatory data. In reference 2, the parameters in (6) were consid- 
ered as functions of the angle of attack for the identification of nonlinear models from wind tunnel 
dynamic data with large amplitudes. In both cases, identified models were found to be adequate repre- 
sentations of experimental data. 


A simple form of the indicial function C L{ {t) defining the downwash effect on this tail can be 
obtained by considering that a sudden change m the angle of attack of the wing results in a sudden 
change in the downwash. Downwash from the wing is delayed from reaching the tail by the time 



where Jl f is the distance from the ax. of the wing to that of the tail and V is the airspeed that remains 
constant and equal to the speed of undisturbed flow. By neglecting the unsteady aerodynamic effect on 
the tail, the indicial function C , f (t) can be expressed as 


C Lt (0 = C Li («>)<* -A/) (7) 

Its time history is sketched in figure la. This assumption about the downwash has been used by Cowley 
and Glauert in formulating contributions of the tail to the acceleration stability derivatives C L and 


A more detailed qualitative assessment of C } (t) is given in reference 6. At the time of sudden 
change in a u , and under condition a ( = 0, no disturbances are present at the tail. At certain time 
> 0 , the tail begins to penetrate the wake created by the vorticity shed by the wing. The tail begins to 
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develop positive lift on that portion of its surface that has penetrated the wake. As time increases the tail 
further penetrates the field of shed vorticity and is subjected to a combined effect of upwash and down- 
wash. Further penetration into the wake results in a change of lift from positive to negative. At certain 
time t 2 >t\ the tail is fully immersed in the steady downwash field. Then the lift on the tail remains 

constant because it is no longer influenced by the disturbances created by the wing at the beginning of 
the motion. 

From the preceding discussion one can anticipate the shape of C u (/) as sketched in figure lb 
and its analytical form as a combination of three effects: 

a) change in the lift of the wing followed by a change in angle of attack of the wing; 

b) change in the downwash at the tail following a change in lift on the wing; 

c) change in the lift on the tail due to a change in downwash at the tail. 

Therefore, the variation in the lift on the tail following a change in angle of attack of the wing, while 
maintaining a t = 0 , follows as a solution of three equations 


C LwW ~ C Lw (0«(°)+ fc LH , {t-x)a{x)dx (8) 

0(H J OCVV 

0 


£(0 = e CL (?)C LM .(0) + Je CL^ r ~ x )^Lw( x )dT (9) 

0 


C Ll (t) = - c u (?)e(0) - f C L (t-x)t{x)dx (10) 

K J # 

0 

The minus signs in equation (10) reflect the definition of a positive downwash angle that results in a 
decrease of the angle of attack at the tail. The corresponding indicial function C. (?) can be obtained 
from equations (8) to ( 1 0) as a " 


£ a (?) " £ CL^ C Lw a J 0 ^ + - x ^ LWn JX)dx (||) 

0 


C uJ» = - c l> (Oe a (0)- \c Ll (t - x)t a (x)dx (12) 

0 

where the last equation can be also written as 
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( 13 ) 


C L , JO = -ea(')C 11|i (0)-Je a (»-T)C t ,,(T)<ft 

0 


In the above equation the indicia! function (?) has already been discussed. The function 

represents the lift on the tail during passage of the tail through a step change in downwash that is equiv- 
alent to the passage through a sharp-edge gust. The two-dimensional theory for this case was developed 
by Kiisner in reference 9 and extended to the three-dimensional case in reference 10. It was found that 
the indicial function C L (?) can again be approximated by a sum of exponentials similar to equation 

(5) for C L (?) . 

Analytical expressions for the indicial function e CL (?) in equations (9) and ( 1 1 ) were developed by 
several authors, see e.g., references 5, 7, and 8 using various degree of approximation to the vortex 
sheet behind the wing. A rather simple approach to this problem is demonstrated in the next section. 
Once the indicial function explaining the variation in lift of an aircraft following changes in the angle of 
attack are known, they can be converted to those in the pitching-moment equation. 


Downwash Angle 

In the following development the vortex sheet of a finite wing in an incompressible steady flow is 
replaced by a horseshoe vortex consisting of the bound vortex and two trailing vortices. After a sudden 
change in the circulation around wing, a stationary vortex separates from the trailing edge of the wing 
and moves in the direction of the velocity vector with the velocity equal to that of free stream, see 
figure 2. The vortex system shown in figure 3 indicates the position of a reference point A and distance 
between trailing vortices. For the present study this distance is equal to the wing span. As pointed out in 
reference 1 1, however, the distance between trailing vortices is less than the wing span. It can be deter- 
mined by equating the lift of the replacement system to the lift of the wing. 

The downwash angle at the point A is £ = w/V where w ( is the sum of induced velocities by all 
four vortices. These velocities are calculated in the plane of the vortex system, midway between the two 
trailing vortices, from the Biot-Savart rule. The expression for the downwash angle takes the form 



where AT is the circulation strength of each four vortices and L = 3/4 c+i. It is convenient to 
express the downwash angle in nondimensional form as a function of time. This can be achieved by 
considering the relation 


x = Vt-1 


(15) 


where J is the distance between the wing trailing edge (the origin of the wing circulation) and the lead- 
ing edge of the tail. By introducing nondimensional time 
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( 16 ) 


/' 



and nondimensional lengths in terms of the half span, i.e.. 


f = J- and L' = 
b/2 b/2 b/2 


equation (15) takes the nondimensional form 


x' = I) 

Substituting (17) into (14) and replacing L by L' yields 


(17) 


£(/') 


Ar 

nVb 


W-i) + 


I 

1) 




J 


1 +L’ 


(18) 


It is further assumed that the Zhukovski’s formula 


pVA Vb = ipK 2 SAC L 

can be also applied to the unsteady flow around the wing. This leads to an increment in the circulation 


Ar 


sv A _ 

2bA ACL 


(19) 


2 

where A = b /S is the wing aspect ratio. Combining (18) and (19) gives the expression for the down- 
wash indicia! function 


E cP'') = 


2nA 


' ru '- ,)+ nh] L ' + h 


i 1 > 2 JT+ 


L’~ 


and it’s steady value 


( 20 ) 




27lA 


I + 


L, + P 


1 n 


J 


l + V 


( 21 ) 


The indicia! function £^(0 is sketched in figure 4. It shows infinite values for t' = 1 which is the 
time when the leading edge of the tail reaches the position of the starting vortex. After that, the down- 
wash angle decreases and, with increasing time, approaches a steady value. Although the indicial func- 
tion can reach an infinite value, the convolution with the lift results in finite downwash at all points. If 
the wing wake passes either above or below the tail surface, the indicial function e c ^(t') will not have 
infinite values. Using the Cowley-Glauert explanation of the downwash delay and considering L = 1 , 
the expression for the downwash indicial function is simplified as 

= ( 22 ) 
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or 


e* CL (0 = £ C l(° o)1 ( r ' _ i) (23) 

where and are the unit step functions. The time history of £* ct (t') is included in 

figure 4 as a dashed line. 

A difference in calculated downwash using indicial function given by equation (20) or (23), and the 
effect of geometry and flight condition will be demonstrated on the harmonic motion of an aircraft. 
Equation (9) has its frequency domain equivalent 

i(ik) = ikt CL (ik)Ciy l /{ik ) (24) 

which leads to the frequency response function 

G Aik) = = ikl CL {ik) (25) 

e C LW (ik) 

where k = coi'/ V is the reduced frequency and 

oo 

Z CL (ik) = je CL (t')e lkt dt' (26) 

0 

For the numerical evaluation of the Fourier integral in (26) it is advantageous to express the indicial 
function e CL (r') in the form 


Z CL (n = Ecz/ 00 )-^'') 

where the deficiency function follows from (20) and (21) as 


F At') = 


2nA 


Ji 


(27) 


Using equations (25) and (27) the frequency response functions relating the downwash and the lift takes 
the form 


G £ (ik) = e CL (oo )-ikjf t it')e ik, 'dt' (28) 

o 

The computed frequency response curve 2nAG e (ik) for £' = 1.0 and L’ = 1.2 is plotted in figure 5 
as a vector diagram. In the same figure the frequency response curve based on the Cowley-Glauert’s 
“lag-in-downwash” formulation 

G*(ik) = e CL (oo)e~‘ kL/I (29) 
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is also given. For the presentation of differences in both downwash formulations and their variation 
with the frequency and distance 1' , the transfer functions G £ (ik) and G* £ (ik) defined by the previous 
two equations are expressed in terms of amplitudes and phase angles as 

G £ (ik) = A t e ,(? e (30) 

G* e (ik) = A\e"*l (31) 

Then the amplitude ratio A £ /A* for three different values of £' = 1.0, 0.5 and 0.25 are plotted in 
figure 6. Figure 7 shows the time delay between harmonic responses obtained from equations (28) and 
(29) as 


<Pe L 


A i 


(32) 


Considering the results in figures 5 to 7, the differences in both downwash formulations depend mainly 
on two parameters, the reduced frequency. A, and nondimensional distance i'.The effect of nondimen- 
sional tail arm V = J!' is rather small as follows from equations (20) and (21). The differences in 
amplitude ratios increase with increased reduced frequency. Increased reduced frequency corresponds 
with rapid maneuvers at low airspeed for an aircraft with long tail arm. On the other hand, these differ- 
ences at given reduced frequency decrease with increased parameter As to the differences in the time 
delays, an increase in reduced frequency makes these differences, in general, smaller. At a given 
reduced frequency the differences in time delays will decrease with the parameter £' which is especially 
apparent at low values of A, see figure 7. More detailed assessment of both downwash formulations can 
be made only for specific aircraft geometry and flight conditions. 


Lift and Pitching Moment 

Returning to equation (4) the lift and pitching-moment coefficient are expressed as 


t 

C L (t) = C L (0) + C L Joo)a(/) + ±;C L </(r)-jF La (r-T)a(T)c/T (33) 

0 

t 

= C m {0) + C WJ ^(°°)a(0 + ^yC m q{t) - jF m Jr-z)aiz)dz (.14) 

0 

Contributions of the wing and tail to these coefficients can be found from the relations 

C L (t) = C Lw (t)- S jc u (t) (35) 

CJt) = C mi + hC L (t)-V l C u (t) (36) 


where C is the pitching moment about the aerodynamic center of the wing, ax., h is the distance 

tit . 

between the ax. and the aircraft center of gravity expressed as a fraction of c , V f = f /Sc is the tail 
volume and C Lt is based on S { . 
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Combing equations (33) through (36) yields 



C L (t) = C L {0) + C L "(oo)a(t) + ^pC L (oo)q(t) 


/ 

-1 

0 

s s 

f l» V~ x ) + T F Lt (f-t) + -r'F L (/-X) d{x)dx 

LW an O S ^avr 

(37) 


+ c , .„(“)«(» + ^ c »„(~>««> 


/ 

-J 

c 

1 

(38) 

where the deficiency functions are defined as 



F L W (0 = C Lw (°°) — ^ Lw (0 

L aw LW aw LM aw 

(39) 


F Lt ^ = C Lt (°°) “ ^ Lt 

(40) 


F Lt = C Lt (°°) “ C \t (0 

^ l aw L 'a« 

(41) 


As mentioned earlier, the steady parts of indicial functions in equations (37) and (38) are aircraft stabil- 
ity derivatives C L , C m , C L and C m . The first two are usually formed as 



(42) 


(43) 


where the contributions of the wing and tail are easily recognized. In these expressions 
a = Ci (oo) and a t = C,, (^o) are the lift curve slopes of the wing and tail respectively. 

H LV au 1 t 

The derivatives C L and C m represent aerodynamic effects that accompany rotation of the aircraft 

about its lateral axis while a remains constant. Both the wing and tail are affected by the rotation but the 
wing contribution is usually negligible. A discussion of the damping term C m due to wing in super- 
sonic flight conditions is contained in reference 6. The contribution of the tail includes changes in the 
lift on the tail due to increased angle of attack by (qJ t )/V and changes due to effective angle of attack 

at the tail caused by the pitching velocity of the wing. The latter change in the lift is discussed in refer- 
ences 12 and 13 for supersonic flight regimes. Thus considering a contribution of the tail only, the 
damping derivatives are approximated as 


C L~Wl 


(44) 


12 


( 45 ) 


Examples 

In the following examples, the indicial functions C L (r) and C m (t) were computed using two sets 
of data representing a fighter and transport aircraft. The geometrical characteristics of these aircraft are 
summarized in table 2. The approximate expressions for the indicial functions of the wing and tail, 
C Lw ar| d C Lt (0 * and the gust lift function C Ll ( t ) were taken from reference 14. The func- 
tionsTor given aspect ratios and argument t' = V t /J ha^e the form 

for fighter aircraft: 

C, (/') = 3.77(1 - 0.283 f -0 ' 626 ') 

U n a vi 

C,. (t r ) = 4.65(1 -0.361 e -0442 ') 

^ at 

C Lt (?') = 4.65(1 - 0.448^’ -<) "^ f,,, - 0.272e _<) ' 841 ^ - 0. 1 93e ) 

for transport aircraft: 


C., ({') = 4.65(1 - 0.361 e Xi2! ' ) 

C,. ((') = 3.77(1 -0.283 c -2 ' 20 ') 

Ul a t 

11/ 1 O C 4 f 

C u (t') = 3.77(1 -0.679c -0.227c ) 

* 

The indicial functions e CL (r') and e CL (t') were computed from equations (20) and (23), indicial func- 
tion C (/') from equations (II) and (13). The resulting indicial functions C, (t') and C (/') 

ran L a m u 

were obtained from equations (35) and (36) as 

C L^) = C L W ( ( ') + % C Lt ( ? ') + C L, (''» (46) 

a O Ul ai L 'aw 

c m (*') = hC L U')-V T (C Ll (t') + C L (/')) ( 47 ) 

a 1 u l at L/ avt 

assuming C m = 0 . 

The three indicial functions of the fighter aircraft appearing on the right side of equations (46) and (47) 
are plotted in figure 8, the functions C L (t’) and C m ( 1 ') in figure 9. The dashed lines in these figures 

indicate the “lag-in-downwash” approximation of the indicial function e CL (t'). Figure 10 and 1 1 

present the same functions as above for the transport aircraft. As expected, the form of the indicial func- 
tions C L (/') for both aircraft was only slightly affected by different expressions for e r/ (/') . This 

a k l 

means that the model structure for C L ( t ') can be postulated as that for the wing alone, i.e.. 
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C L a {t) = C o + 1 J C / J ' 

The structure for an adequate model for C m (r'), on the other hand, could be strongly affected by 
t CL (t'), evaluated either from equation (20) or (23). As mentioned previously, the decision between 
using one of the forms for t CL (t') will depend on the values of two parameters, k and A' . Assuming a 
harmonic motion with co = 4 rad/sec and V = 1 00 m/sec, and given geometry, then for the fighter aircraft 
k = 0.08 and A' = 0.44 . Using figures 6 and 7 it can be found that the amplitude ratio, A £ ~ 1.0 and 
time delay 5 = 0.01 7 sec. In this case, therefore, the Cowley-Glauert approximation for e CL would 
be warranted. For the transport aircraft k = 0.96 and A' = 0.80 which resulted in A £ ~ 2.0 and 
5 = 0.12 sec. In this case, equation (20) for e CL (r') should be used. Then the model for C m (/')canbe 
postulated as 

£„„(') - ‘V rf i<'-' 1 >“ < 49 > 

or 

C m U) = d 0 + d + d 2 (t - 1 ] ) + + d^(t - 1 2 ) + (50) 

The plus functions in (49) and (50) are defined as 

0 lfor 'l^' 

(t-t,) = 

+ 0 tor 1 1 < t 

t - 1 j for / j > t 


and similary for (r - r 2 ) + • 

Equation (49) reflects sudden change in the downwash angle, whereas equation (50) takes into account 
gradual changes in the downwash. In equation (49) the value of d 0 would be different from that based 
on Cowley-Glauert’s approximation of z CL (t) ■ The time r, in both expressions might be close to 
t = A/V. 

Concluding Remarks 

Aerodynamic model equations for the longitudinal motion of an aircraft with horizontal tail were devel- 
oped. The equations include unsteady effects modeled by linear indicial functions. These functions rep- 
resent responses in the lift and tail alone, and interference between those two lifting surfaces. The 
indicial functions for the lift on the wing and tail were modeled by a series of exponential functions. The 
effect of the wing on the tail was formulated for two different expressions for the downwash angle at the 
tail. The first expression used the Cowley-Glauert approximation known as “lag-in-downwash.” The 
second took into account growth of the wing circulation and delay in the development of the lift by the 
tail. From analytical expressions obtained and numerical examples, the following conclusions can be 
drawn: 

a) for harmonic motion the differences in calculated downwash using simple “lag-in-downwash” 
concept and that with unsteady aerodynamics depend upon the frequency and airspeed of the 
motion, on the distance between the wing trailing edge and tail leading edge, and on the wing 
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span. These differences will increase for an aircraft with long tail arm performing rapid maneu- 
vers at low speed. More detailed assessment of both concepts, however, can be made for speci- 
fied aircraft geometry and flight conditions; 

b) a model for the lift is only slightly affected by different approaches in computing the downwash 
angle. For the preliminary studies and model identification, the model structure can be the same 
as that for the wing alone; 

c) in pitching-moment formulation the “lag-in-downwash” approximation is warranted for aircraft 
with short tail arms. For aircraft with long tail arms, however, the expression with unsteady 
terms should be considered. Then the resulting indicia! function for the pitching moment can be 
approximated by zero-order or first-order splines. 
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Table 1 . Summary of components of indicial functions expressing changes in lift due 

to unit step in angle of attack. 



Table 2. Geometrical characteristics of two 
aircraft. 












c. 


C, 



(b) Gradual vortex sheet penetration. 

Figure I . Qualitative time histories of tail lift caused by unit step in wing angle of attack 



Figure 2. Location of bound and starting vortices. 
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Figure 5. Frequency 
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Figure 8. Indicia] functions of wing and tail. Fighter aircraft. 
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